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Abstract 

In cosmological first-order phase transitions, the microscopic interaction of the 
phase transition fronts with non-equilibrium plasma particles manifests itself macro- 
scopically as friction forces. In general, it is a nontrivial problem to compute these 
forces, and only two limits have been studied, namely, that of very slow walls and, 
more recently, ultra-relativistic walls which run away. In this paper we consider 
ultra-relativistic velocities and show that stationary solutions still exist when the 
parameters allow the existence of runaway walls. Hence, we discuss the necessary 
and sufficient conditions for the fronts to actually run away. We also propose a phe- 
nomenological model for the friction, which interpolates between the non-relativistic 
and ultra-relativistic values. Thus, the friction depends on two friction coefficients 
which can be calculated for specific models. We then study the velocity of phase 
transition fronts as a function of the friction parameters, the thermodynamic pa- 
rameters, and the amount of supercooling. 

1 Introduction 

Cosmological phase transitions may have observable consequences. In particular, first- 
order phase transitions provide a departure from thermal equilibrium, which may give 
rise to a variety of cosmological relics, such as the baryon asymmetry of the universe 
[1], cosmic magnetic fields |2], topological defects |3], baryon inhomogeneities |11[5], and 
gravitational waves [6]. In a first-order phase transition, bubbles of the stable phase 
nucleate and grow inside the supercooled phase. The velocity of the phase transition 
fronts (bubble walls) is an important parameter for the generation of cosmological relics. 
For instance, generating sizeable gravitational waves requires high velocities, whereas the 
generation of baryon number in the electroweak phase transition peaks at small velocities. 

The wall velocity is governed in principle by the pressure difference between the two 
phases. However, the wall propagation disturbs the plasma, and the latter resists the 
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motion of the wall. This opposition manifests itself in two ways. Microscopically, the 
interactions of non-equilibrium plasma particles with the wall cause a friction force on the 
latter [U El El [10] • Besides, the latent heat that is released at the phase transition fronts 
causes temperature variations in the plasma [HI |T2l [131 [E] ? which tend to diminish the 
pressure difference between phases. These two mechanisms are called microphysics and 
hydrodynamics, respectively. 

Microphysics is a very difficult subject and for several years was considered only in 
the non-relativistic (NR) limit. As a consequence, phenomenological models have been 
used in order to extrapolate the friction force to higher velocities. The general approach 
is based on adding a covariant damping term rju^df^cj) to the field equation, where u'^ 
is the four- velocity of the plasma (see, e.g., [5l [12]). This gives a friction force of the 
form Ffr ~ rj'y^Vw, where is the wall velocity and 7^ = 1/ ^/T^^v^. The parameter 
7] can be determined by considering the limit f ^ — )■ and comparing with microphysics 
calculations. 

Recently, the opposite limit was considered [I0|. The total force acting on the elec- 
troweak bubble wall was derived for an ultra-relativistic (UR) wall propagating with 
extremely large values of the gamma factor. For such a fast moving wall the physics is 
quite simpler than in the NR case, since the plasma is almost unaffected by the passage 
of the wall. The resulting total force does not depend on the wall velocity. As a conse- 
quence, if the force is positive, then the wall runs away (i.e., as the wall propagates, 7^ 
grows linearly with the propagation distance). If the total force turns out to be negative, 
then the initial assumption of an extremely ultra-relativistic wall is incorrect and the wall 
must reach a slower, stationary state. Thus, we have a simple criterion for the existence 
of runaway solutions. The results of Ref. [TU] motivated modifications to the usual phe- 
nomenological models for the friction, so that the friction saturates for large 7^ [T5l [T6] . 
The modifications essentially amount to eliminating the gamma factor from the friction 
force, so that Ffr ~ r/f^. 

We wish to point out that the runaway condition found in Ref. [lOj is only a necessary 
condition for the existence of the runaway solution. It does not guarantee that the runaway 
solution will actually be realized. Indeed, it is well known that different propagation modes 
can exist for the same set of parameters (due to nonlinear hydrodynamics). In particular, 
the runaway solution may coexist with a detonation solution. Since detonations are 
presumably stable, the wall probably will not run away whenever it can propagate as a 
detonation. We also notice that a single-parameter phenomenological model of the form 
Ffr = rjVw will hardly give the correct values of the friction force in both ultra-relativistic 
and non-relativistic limits. 

The main goals of this paper are the following. In the first place, we wish to identify 
the friction force in the UR limit, in order to treat very fast but stationary solutions. This 
can be done by decomposing the total force obtained in Ref. pTO] into a driving force and a 
friction force. In principle, though, it is not clear which is the correct decomposition, since 
there is no velocity-dependent term in this force. Besides, hydrodynamics is different in 
the stationary and runaway regimes. We wish to consider the existence of stationary and 
runaway solutions and discuss the necessary and sufficient conditions for the wall to run 
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away. Secondly, we wish to construct a well-motivated phenomenological model, which 
interpolates between the NR and UR limits. The model will have two free parameters 
and will allow us to treat stationary solutions with intermediate velocities. We also wish 
to discuss the case of phase transitions with a large amount of supercooling, which favor 
fast-moving phase transition fronts. 

The plan is the following. In section |2] we briefly review the dynamics of cosmological 
phase transitions and the hydrodynamic solutions for stationary fronts. In section [3] we 
discuss the friction forces. We introduce an ultra-relativistic friction parameter and we 
discuss the conditions for the wall to run away. In section S] we discuss a phenomenological 
model for the friction. In section [5] we consider the bag equation of state to obtain 
analytical formulas for the wall velocity and for the runaway conditions. We also analyze 
the general dependence of the bubble wall velocity on the friction and thermodynamic 
parameters. Finally, in section E] we discuss the case of strong supercooling. We summarize 
our results in section [71 

2 Phase transition dynamics 

A cosmological phase transition is described in general by a scalar field which acts 
as an order parameter. We shall consider the case in which is a Higgs field. The 
finite temperature behavior of the system is determined by the free energy density (finite- 
temperature effective potential) 

^(0,T) = V(0) + VV(0), (1) 

where V{(j)) is the zero-temperature effective potential and Vt{4>) the finite-temperature 
correction. To one-loop order, the latter is given by [T7] 

VTi<P) = Y:±9.f^,lo,{lTe--^n- (2) 

where the sum runs over particle species, Qi is the number of degrees of freedom of 
species i, the upper sign stands for bosons, the lower sign stands for fermions, and 
Ei = a/p^ -l- mf{(j)). Here, rrii are the field-dependent masses. There is also a correction 
from the resummed daisy diagrams for bosons. For our discussions it is not necessary to 
consider the exact form of this correction. Its effect is to modify the cubic term in the 
expansion of Vr in powers of rrii/T, as we shall comment below. 

In thermal equilibrium, the field lies at a minimum of the free energy. If there are 
several minima, then different phases are possible for the system, and phase transitions 
may occur. In the simplest case we have a high-temperature minimum 0+(T) (in general, 
0+ = 0) and a low-temperature one 0_(T). For a first-order phase transition, these 
two minima coexist in a certain range of temperatures, separated by a barrier. All the 
properties of a given phase are derived from J-", once it is evaluated at a minimum. Thus, 
the two phases are characterized by the free energy densities 

J-+ (T) ^ ^(0+ (T) , T) , (T) ^ ^(0_ (T) , T) , (3) 
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which give different equations of state (EOS). The energy density is given by p± (T) = 
J-'±{T) — TJ''_i_(T), where a prime indicates a derivative with respect to T. The pressure 
is given by p±(T) = —J^±{T). The enthalpy density is given by w± = p± + p±, and the 
entropy density by s± = w±/T. The speed of sound is given by c^(T) = dp±/dp± = 
p'-i-{T) / p'_^(T) . The phases are in equihbrium at the critical temperature Tc, defined by 
^+(Tc) = J-'_(Tc). The energy density difference at T = Tc is called the latent heat 
L^p+ (T,)-p_ (T,). 

The phase transition occurs, in principle, when the Universe reaches the critical tem- 
perature. However, the nucleation rate vanishes at T = Tc. Bubbles of the stable phase 
effectively begin to nucleate at a temperature T„ below Tc [18j. The nucleated bubbles 
expand due to the pressure difference between the two phases (at T < Tc). The phase 
transition fronts move with a velocity v^j which depends on the pressure difference and 
the friction with the surrounding plasma. The release of latent heat causes local reheat- 
ing and bulk motions of the plasma. Considering the development of a phase transition 
involves solving a set of coupled equations for several quantities, each of which is not easy 
to compute, such as the nucleation rate, the wall velocity, and the temperature. The 
latter varies due to the expansion of the universe and the release of energy at the phase 
transition (for the dynamics of phase transitions see, e.g., [12]). Here we shall concentrate 
on the motion of phase transition fronts for a given nucleation temperature T„. 

Thus, the relevant variables are the scalar field 0(x, t), the temperature T(x, t), and 
the velocity v (x, t) of the plasma. We shall consider the thin- wall approximation for the 
field profile. Thus, we shall assume that the field varies in an infinitely thin region, outside 
which is a constant. For the macroscopic treatment this is a good approximation, since 
the wall width is much smaller than the width of the fluid profiles [20] • We shall consider 
planar-symmetry fronts moving in the z direction (see [21] for a discussion on considering 
different wall geometries). As a consequence of the friction with the plasma, the bubble 
wall often reaches a terminal velocity in a very short time after bubble nucleation. In 
the rest of this section we shall review the hydrodynamics for stationary phase transition 
fronts. 

2.1 Hydrodynamics for stationary fronts 

The equations for the fluid variables can be obtained from the conservation of the energy- 
momentum tensor, d^T'^'^ = 0. For stationary profiles moving at constant velocity, it is 
useful to consider the rest frame of the front. The relation between the fluid variables on 
each side of the wall is well known |22j , 

t(;_t>_7^ = w+i;+7^, (4) 
w_vt'^t+P- = w+f^7^+p+, (5) 

where 7± = 1/ a/1 — f^, and we have used a -|- sign for variables in front of the wall 
and a — sign for variables behind the wall (which correspond to the -|- and — phases, 
respectively). For an infinitely thin interface, Eqs. (jUlSj) determine the discontinuity of 
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the fluid profiles at the phase transition front. There can also be discontinuities in the 
fluid profiles away from the bubble wall. These are called shock fronts. In such surfaces, 
the temperature and fluid velocity are discontinuous but the EOS is the same on both 
sides. In the reference frame of shock fronts Eqs. (jMS]) still apply. 

Since there is no characteristic distance scale in the fluid equations, it is usual to 
assume the similarity condition [22], namely, that w,p and v depend only on ,^ = z/t. 
In the planar-symmetry case we have either constant solutions v{^) = constant, or a 
"rarefaction wave" solution frar(0 = ~ c)/(l — C,c), where c is the speed of sound (see, 
e.g., [211 120] for recent discussions). The fluid velocity profile is constructed from these 
solutions, using the matching conditions (jUlS]) and appropriate boundary conditions. 

For a given set of thermodynamical parameters, Eqs. (jl]) and (|5]) give f + as a function 
of u_. The solutions have two branches, called detonations and deflagrations (see, e.g., 
[23] )• For detonations the incoming flow is faster than the outgoing flow (|f+| > l^-l) and 
is supersonic (|f+| > c+) for any value of f_ in the range —l<v^< 0. For deflagrations, 
we have If+I < and |f+| < c^. At = c_, the value of |f+| is a minimum for 
detonations and a maximum for deflagrations. At this point the hydrodynamic process 
is called a Jouguet detonation or deflagration. We denote the corresponding values by 
|f+| = ^5^* and If+I = f^'^^, respectively. The hydrodynamic process is called weak if 
the velocities v+ and f_ are either both supersonic or both subsonic. Otherwise, the 
hydrodynamic process is called strong. 

The fluid temperature and velocity profiles depend on the boundary conditions far 
behind the wall (at the center of the bubble) and far in front of the wall, where information 
on the bubble has not arrived yet. Thus, the fluid velocity vanishes far behind and far in 
front of the wall. Furthermore, the temperature far in front of the wall is given by the 
temperature T„ at which the bubble nucleated. Three kinds of solutions are compatible 
with these requirements: a weak detonation, for which the wall is supersonic, a Jouguet 
deflagration, which is also supersonic, and a weak deflagration, which is subsonic. 

The relevant properties of these solutions are the following (see, e.g., [201 EI] for 
details). For the supersonic detonation, the fluid in front of the wall is unperturbed and 
we have Vw = — t'+,T+ = T„. Hence, we have Vw > vf^^. The fluid profile behind the 
wall is given by the rarefaction wave. The strong detonation solution must be discarded 
because it is not possible to construct a fluid profile for it. For the deflagration, the fluid 
velocity in front of the wall is not at rest, and there is a supersonic shock front preceding 
the wall. Beyond the shock, the fluid is still unperturbed at temperature T„. The solution 
between the phase-transition and shock fronts is a constant. The boundary condition for 
deflagrations is often considered to be that the fluid behind the wall is at rest. In this case, 
we have Vu, = —V-. This "traditional" deflagration can be either weak (and subsonic), or 
strong (and supersonic). The limit between these solutions is a Jouguet deflagration with 

= c_. It has been argued that strong deflagrations are unstabl^ and should not be 
considered [121 [131 E] • However, supersonic deflagrations that are not strong can exist 
[13] if we relax the boundary condition that the fluid is at rest behind the wall. Thus, if 

^Even weak deflagrations may be unstable 
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the condition = —v^ is replaced by f„ = — c_, we have a Jouguet deflagration with 
= Vj^^ < c+. In this case the wall is followed by a rarefaction wave. Since the wall 
moves at the speed of sound with respect to the fluid behind it, and the fluid also moves 
with respect to the center of the bubble, the wall velocity is supersonic. This solution fills 
the velocity gap between the weak deflagration and the detonation, C- < < Vj^^. 



3 Microphysics 

The forces acting on the bubble wall can be derived from the equation of motion for the 
background field 0. It is usual to consider the WKB approximation, which makes sense 
if the scale of variation of is not too short. We have [7j 

where fi is the distribution function of particle species i. In general, a distribution function 
can be written as 

m = frip)+5Mp,x), (7) 

where fi'^{p) = l/(e^*/"^=Fl) is the equilibrium distribution function and Sfi is a deviation. 
Inserting this decomposition into Eq. (|6]), the term fP just gives dYr/dcl) (i.e., the finite- 
temperature corrections to the effective potential) and we obtain 

We shall assume for simplicity that the phase transition occurs from a vanishing minimum 
to a non- vanishing one, as in the usual symmetry braking case. We shall denote 0o the 
non- vanishing minimum, i.e., 

0+(r)^O, <P4T)^MT) (9) 

An equation for the bubble wall can be obtained if we assume a fixed field profile 
centered at the wall position Zu,(t). Thus, consider the dependence (j){z, t) = (f)['jy^{z — Zu,)]. 
The exact form of the one- variable function (j){z) is not relevant. It is often assumed to be 
given by a tanh. It varies from the value 0(— oo) = 0o (inside the bubble) to 0(+oo) = 
(outside the bubble). The variation occurs in a small region of width (the bubble wall). 
In the reference frame which instantaneously moves with the wall, in which z^, = and 
z^ = 0, we have 7^ = 1,'jw = and = where a dot indicates a derivative with 
respect to t. The first term in Eq. (|8]) gives [z^^z — Zy,)(f)'{z) — <f)"{z), where a prime 
indicates a total derivative with respect to z. Notice that the function (^'{z) has a peak 
inside the wall and vanishes outside. We can define the wall position so that J z(f)''^dz = 0. 
Then, if we multiply Eq. by (p' and integrate across the wall, we obtain 

az^ = FdJA + FiJA, (10) 
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where o' = J (j)''^dz is the surface tension, and 

A J d(j) dz 



dz, 



(11) 



^fr / , am- 

T = z:^. dz^ 




(12) 



are the forces per unit area acting on the wall. 

The force Ff,. depends on the interactions of the plasma particles with the wall and may 
be regarded as a friction force. Indeed, for small wall velocity (in the reference frame 
of the bubble center), this force turns out to be proportional to Vw As a consequence, 
the wall may reach a terminal velocity. The force does not depend on microphysics 
details, and may be regarded as the driving force. 

3.1 The driving force 

To see the behavior of this force, let us consider a wall which moves very slowly, so that 
equilibrium can be assumed and Ffr vanishes in Eq. (ITU]) . For such a slow wall, the 
temperature will be homogeneous, T{z) = constant. Hence, the integral in Eq. fllip 
yields J^{(t>+, T) — J^(0_, T). Thus, we obtain the equation = P-{T) — p+(T), which is 
positive for T < Tc. Therefore, if the wall is initially moving very slowly, it will accelerate 
due to the pressure difference. As a consequence, the wall motion will cause departures 
from local equilibrium and a friction force will appear. Besides, inhomogeneous reheating 
will arise, which affects the driving force. 

For the integral in Eq. flTT]) we may use the identity (dJ^/d(l)){d(l)/dz) = dF/dz — 
{dJ^/dT){dT/dz) |l2l|26], so that the driving force can be expressed as 



where the entropy density is given by s = —dT/dT. In order to obtain analytical results, 
in this paper we shall use approximations for integrals across the wall. Our aim is to 
obtain expressions which only depend on variables defined outside the wall, namely, the 
values of v±, T±, etc., which can be obtained using Eqs. (HHS]), the EOS, and appropriate 
boundary conditions. For the integral in Eq. f|T3|) . we shall approximate the entropy 
density by the average value (s) = (s+ + s_)/2. This gives the approximation [26j 



If we neglect hydrodynamics and consider a homogeneous temperature, we obtain the 
pressure difference P-{T) — p^{T). This is very sensitive to the departure of T from T^. 
Besides, p^(T_) — p+(T+) is very sensitive to the difference T+ — T_. Indeed, the pressure 
difference may be positive or negative depending on the value of — T_. The general 
effect of hydrodynamics is to slow down the wall [26| |27] . 




(13) 



p_(r_)-p+(r+) + (.) (T+-r_). 



(14) 
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3.2 Friction force 



The friction force Ff^ is in general much more difficult to calculate than F^,.. It has 
been extensively studied for small wall velocities [71 18] , since temperature gradients can 
be neglected and the deviations from equilibrium can be assumed to be small. Even in 
this case, computing the deviations 6fi involves solving a complex system of Boltzmann 
equations. Recently [lU], the total force acting on the bubble wall was derived for a wall 
which propagates ultra-relativistically. The treatment turns out to be quite simpler. 



3.2.1 Non-relativistic limit 

For small velocities we have 6fi oc Vw, and Eq. (fT2|) gives a friction term of the form 

^ = -VnrVw (15) 

The non-relativistic friction coefficient r^NR depends on the couplings of the particles to 
the Higgs, and also on the interactions of plasma particles, which tend to restore the 
equilibrium. Effective interaction rates Tij ~ T appear in the equations for 6fi, coming 
from the collision integral in the Boltzmann equation. The "thick wall" limit F ^ 1//^ is 
often used, which leads to analytical results (see, e.g., [S]). Here we shall use the results 
from the simplified treatment of Refs. [25] . For masses of the form = hfcfP' + /i^ 
we have 

l4 r+oo 

^NR^En^/ ^l-X<P)<P'rdz. (16) 
• J- J-oo 

where F is an average interaction rate (typically, F < lO^^T), and the function ci(0) is 
given by 

1 f <fp e^/'^ 
""'^T^J {27rfE{e^/^Tl?' ^^^^ 

Notice the strong dependence of the coefficient ci on the ratio m/T. Different limiting 
cases have been analyzed in Ref. [28]. For a "typical" phase transition with (p T ^ Tc, 
one can often use an expansion in powers of m/T. To lowest order, we have ci = log x/27r^, 
where x = 2 for fermions and x = T/m for bosons. Hence, in this case the dependence 
on is at most logarithmic, and we can regard constant for the integral in (1161) . 

We obtain 

r^NR = ^ 1 <P'<l>"dz, (18) 

where 

^NR = Ef7T(,^J (19) 

The dimensionless coefficient i7nr depends on details of the model, whereas the inte- 
gral in Eq. f|T8l) depends only on the wall shape. This integral can be estimated as 
J (p'^cp'^dz ~ 0Qcr ~ (po/lw, where Iw is the wall width. Typically, > T~^. These rough 
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approximations give in principle the correct parametric behavior, but omit numerical 
factors which may be as high as ~ 10. To obtain quantitatively useful values, we may 
"calibrate" the friction with a known model [TBIES]- In our approximation (1191) . this can 
be accomplished by choosing a suitable value of F/T. For instance, considering the top, 
Z and W contributions to f|T9l) . we obtain the correct SM values (?7nr ~ 0.6 [16]) if we 
use F ~ lO-^T. 



3.2.2 Ultra-relativistic limit 

For a wall which has reached ultra-relativistic velocities with very large gamma factor, 
several approximations are justified [lOj. In the frame of the wall, plasma particles always 
have enough energy to surpass the wall, and the reflection coefficients are exponentially 
suppressed. Incoming particles have received no signal that the wall is approaching and 
are in equilibrium. Besides, interactions between plasma particles are time delayed and 
can be neglected. Thus, the occupancies evolve undisturbed. As a consequence, only the 
equilibrium occupancies of the symmetric phase are needed in the calculation. To lowest 
order in 1/7^, the net force per unit area acting on the wall is given by [10] 

^ = V{ct>^) - - 5^^7.K'(0-) - ^^{^^)] I J^^^^f^W- (20) 

This force does not depend on the wall velocity. Therefore, if it is positive the wall will 
accelerate indefinitely. On the other hand, if F/A is negative, then the wall in fact cannot 
reach the ultra-relativistic regime. 

Notice that the result f l20|) can be obtained directly from Eq. ([6]), by evaluating 
the momentum integral in the + phase. Thus, it does not require the decomposition 
into equilibrium occupancies and deviations. We are interested in such a decomposition, 
though. More precisely, we intend to decompose the total force into a driving force and a 
friction force. This will help us to introduce, in the next section, a phenomenological model 
for the friction which interpolates between the non-relativistic and the ultra-relativistic 
limits. Adding and subtracting in Eq. fl20p the finite-temperature correction Vt+(0+) for 
the -|- phase, we have 

^ = ^(0+,r+)-^(0_,r+), (21) 

where J-'(0_,T+) is the mean field effective potential, obtained by keeping only the 
quadratic terms in a Taylor expansion of Vr about the + phase [IHl fT5] . 



H<P-,W = V{<P-) + %(0+) + 5^K(0-) - m^(0+)] ^ 



(22) 



In contrast to the non-relativistic case, the total force in (!2T|) does not have a velocity- 
dependent term, which could be identified as a friction force. Nevertheless, we can still 
decompose the force into a part which comes from the equilibrium distributions (the 
driving force) and a part coming from the deviations (the friction force). The particles in 
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front of the wall are in equilibrium at temperature T+ and pressure p+(T+) = — T+). 
The departures from equilibrium occur at the wall and behind. The corresponding friction 
force must come from the term — T_|_) in Eq. f l2ip . We only need to isolate the 
equilibrium part. 

The problem with such a decomposition is that it is not clear which would be the 
temperature behind the wall. In the stationary case the temperature T_ can be calculated 
from Eqs. (jHE]) and is different from T+. In the limit — )■ 1 (i.e., for a very fast 
detonation) we still obtain T_ 7^ T+. Nevertheless, there is no reason to assume that the 
stationary solution should match the runaway solution. Notice that the only temperature 
appearing in Eqs. fl21ll22p is T+, suggesting that, in the runaway case, we have T_ = r+. 
Indeed, for stationary solutions the released latent heat goes into bulk motions of the fluid 
and reheating of the plasma, whereas for runaway solutions the energy goes mainly into 
accelerating the wall. Assuming T_ = T^, the occupancies in each phase (in the plasma 
frame) are given by 

= l/[exp(E,+/r+)Tl], (23) 
f,.{p.) = l/[exp(E,_/T+)Tl] + 5/„ (24) 

with Ef^ = p± + m^±- 

In fact, we can calculate the exact form of the occupancies just behind the wall. After 
the passage of the wall the occupancies fi+{p+) are undisturbed, but the energy and 
momentum of a particle have changed. In the plasma frame, these changes are given by 

m 



1 1 

mf — mf 



We thus have £'j+(p+) = Ei_{p^) + (m^_^ — m^„)/[2(£'i_ — j9_)], and 



2 ^2 \ H -1 



^i-(p-) , mf+ -m,_ 



(26) 



It seems natural to decompose the occupancies fl2Bl) in the form fl2^ . which can be taken 
as the definition of 5fi. If, e.g., the mass difference mj+ — rrii^ is small, then the deviation 
Sfi will be small too. 

According to Eqs. (124|) and (|26|) . the deviations from equilibrium vanish only in front 
of the bubble wall, in contrast with the non-relativistic case, in which the deviations vanish 
also behind the wall. Indeed, the (0')^ factor in Eq. fITB]) indicates that the deviations 
are localized at the bubble wall. The difference arises because, in the UR case, the wall 
has passed so quickly that there was no time for the plasma to recover the equilibrium. 
Notice, anyway, that Eq. fl26|) gives the occupancies behind the wall but close to it. At 
some distance behind the wall the plasma will reach the equilibrium (and some reheating 
will occur). 

It is easy to obtain the macroscopic version of the decomposition ([MD, directly from 
Eq. (12T|) . For r_ = T+, the driving force is given by 

^ = ^(0+,T+)-^(0_,T+). (27) 
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Therefore, the departure from equihbrium causes a friction force given by 

^ = ^(0_,T+)-^(0_,T+), (28) 

Since ^ 1, we can define an ultra-relativistic friction coefficient t^ur by 

VVK = -FiJA, (29) 

so that the total force is given by 

F 

_ = (p_ - p^) - r/uRf (30) 

Notice that we may also write Eq. (12T]) as F/A = J^(0+,T+) — J^(0_,T+), since 
J-'(0+,T+) = J-'(0+,T+). If this mean-field potential difference is negative, then the 
bubble wall cannot run away. For instance, the wall never runs away in a "fluctuation 
induced" flrst-order phase transition [TU], i.e., a phase transition which is flrst-order due 
to the thermal part of the effective potential, Vr- For example, consider the well known 
high-temperature expansion of Eq. (|2]), 

VV(0) = _^___ + ^^,c,^^ 2^ 9^^^ + Oim), (31) 

i i bosons 

which is valid for many models. Here, Q = 1 (1/2) and = 1 (7/8) for bosons (fermions). 
Only bosons contribute to the rnf terms. In fact, after the resummation of daisy diagrams, 
only the transverse polarizations of gauge bosons remain in general in this contribution. 
The cubic term in Eq. f l?!]) is very important since it may cause a barrier in the free 
energy and, thus, a flrst-order phase transition. This term is not present in the mean fleld 
potential. If the flrst-order character of the phase transition is due to this term alone, 
then in the mean fleld potential the minimum 0_ will raise above the minimum 0+. As a 
consequence, the net force fl2ip will be negative and the wall will reach a terminal velocity. 

Let us consider the friction coefficient ?7ur for the case m/T <^ 1 (as we did for ?7nr)- 
According to Eq. f l3T|) we have, to lowest order in m/T, 



VVR 



bosons 



E f§K(<^-)-^'(<^+)]- (32) 



We have assumed that the contribution of gauge bosons is important. Otherwise, the 
(9(m^) terms should be considered. The generalization is straightforward. Let us further 
simplify the problem by specifying to the case 0+ = and assuming particle masses of 
the form mj(0) = hi(f). Then, we have 

^ E fl^^^^o = VvkT<PI, (33) 

bosons 

where we used again the notation 0o for the symmetry-braking minimum, and we have 
defined a dimensionless coefficient i7ur which contains the parameters of the model. 
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3.3 Runaway conditions 

According to Eqs. f llOp and (1301) . in the runaway regime we have 

az^ = P-{T+) - p+{T+) - r/uR, (34) 
which gives the necessary condition for the wall to run away, 

p_(T+)-p+(T+)>r7uR. (35) 

This condition does not depend on the non-relativistic friction parameter ?7nr and is 
equivalent to J^(0+,T+) — J^(0_,T+) > 0, which is the criterion provided in Ref. [lU]. 
We wish to emphasize that Eq. fl5^ is just a necessary condition: if it is not fulfilled, 
then the wall cannot run away (since the total force in the runaway regime would then 
be negative). However, it is not a sufficient condition for the wall to run away, since it 
was obtained by assuming that the wall is already in the runaway regime. Thus, Eq. (135|1 
implies that the runaway solution exists, but does not guarantee that it will be realized. 
It is well known that multiple hydrodynamical solutions may exist for the same set of 
parameters. Before the wall reaches ultra-relativistic velocities, the forces are not like in 
Eq. f lM|) . and it may be possible for the friction to compensate the driving force. Even 
if the friction force is smaller for smaller velocities, hydrodynamic effects slow down the 
wall, acting effectively as a friction |26] . 

To be specific, consider a very fast detonation. In the limit in which the terminal 
velocity — )■ 1, we have f± — 1. However, for a detonation, the difference of temperature 
between both sides of the wall does not disappear in the limit f„, — > 1. The friction force is 
not too sensitive to temperature gradients, though, and it should be given, for ^ 1, by 
Ff-c/A = T^uR- Thus, the friction for the detonation is as in Eq. (134|) . However, the driving 
force is not, since it is very sensitive to temperature gradients. Using the approximation 
(in]) , we have F^r/A = p_{TJ) — p^(T^) — (s)(T_ — T+). For a stationary solution, the 
driving force equals the friction force, and we have F^r/A = r/uR- For smaller velocities, 
the friction is in principle smaller, and so is the driving force. Hence, for detonations, 
Fdr/A = ?7uR is the maximum value the driving force can reach. We thus have a necessary 
condition for the stationary solution, namely, F^^-fA < 77UR. If this condition is not 
satisfied, i.e., if 

p_(r_) -p+(T+) - (s)(r„ -T+) > r/uR, (36) 

then the stationary detonation solution does not exist. 

We notice that the two conditions ([35D and fl5B]) are really differentia since, for the deto- 
nation, the temperature T_ is always higher than T+. If Eq. ( 136|1 is fulfilled, then the nec- 
essary condition (135|) is also fulfilled. Indeed, given ?7ur and T+, the force p_(T+) —p^(T^) 
is larger than p_(T_) — p+(T+) — (s)(T_ — T+), since the latter is affected by hydrody- 
namics. Therefore, Eq. (1361) gives a sufficient condition for runaway walls, corresponding 
to the detonation velocity becoming v^, = 1. 

^According to the exact expression (fT5)) . for T_ 7^ we could still have F^^/A — p_(T+) — p+(r+) 
if (f){z) = across the wall, which is not the case. 
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According to the above, as the parameters are varied there will exist a runaway solution 
before the stationary solution ceases to exist, and vice versa. Hence, the runaway and 
the detonation solutions will coexist in a certain range of parameters. This range is 
delimited by the conditions fl35|) and f l36|) . Stability analysis indicate that the detonation 
is generally a stable solution [21]. Therefore, we expect that the wall will not runaway 
in the coexistence range. This range will be short if the latent heat is relatively small 
(L -C T^). Indeed, in such a case we will have T+ — T_, and Eqs. and f l5B]) will give 
a similar condition. 



4 Phenomenological model for the friction 

In order to simplify the treatment of phase transition fronts, it is usual to replace the last 
term in the field equation ([H]), corresponding to the deviations from equilibrium, with a 
phenomenological damping term. 

4.1 Existing models 

4.1.1 "Old" phenomenological model 

A widely used approximation is a damping term proportional to u'^d^(j). The coefficient 
of this term can be determined by comparing the resulting friction force (or wall velocity) 
with the results of microphysics calculations such as those considered in the previous 
section. In the general case, the coefficient is field-dependent, 

d,dy+^ + fi(l>)u'^d,(l) = 0. (37) 

We shall refer to this approach as the "old" phenomenological model. 

Proceeding like in the previous section, i.e., assuming a field profile of the form (j){z, t) = 
0[7«)(-2 ~ ^w)], considering a reference frame which instantaneously moves with the wall, 
multiplying by (p' and integrating across the wall, we obtain 

CTZ^ = I ^^dz + j ^vfim'fdz. (38) 



where v{z) is the (negative) fiuid velocity and 7 = 1/a/1 — v"^- The first term in the rhs is 
the driving force, and the last term gives a friction force which, for small velocities, can 
reproduce the form (fT6i) . Indeed, for small the temperature and fluid velocity are not 
significantly altered by the wall. Hence, the temperature is homogeneous and the fluid 
velocity is given by ~ — We thus have 

azy, = P-{T) - p+(T) - r^NRf^,. (39) 

with 

mK = I f{<P)<P"dz (40) 
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The function can be chosen so that Eq. (HOl) reproduces the friction coefficient 

obtained from a microphysics calculation, such as Eq. (1161) . For the case of small rrii/T, 
we can choose 

= Vnr(I>VT, (41) 

so that we obtain the result of Eq. (1181) . 

In the general case, hydrodynamics must be taken into account. In particular, Eqs. 
(jUin]) imply that the fluid velocity and the temperature are in general different on each 
side of the wall. In sec. 13.11 we found an approximation for the first integral in Eq. (1381) . 
For the second integral, we notice that 4''{z)'^ peaks inside the wall and vanishes outside. 
According to the weighted mean- value theorem for integrals [30] , we can replace the func- 
tion f((f))^v with its value at a certain point z = z inside the wall. We shall approximate 
this value by /(0o/2)(7f ), with {jv) = {j-V- + 'y+v+)/2. With this approximation, the 
friction force per unit area, including hydrodynamics effects, is given by 

^ = ^nr(7^). (42) 

For the case of Eq. (HT]) . we have ?7nr = 57nr(0o/4) cr/T. 

Although the phenomenological model (1371) gives the correct form of the friction in 
the non-relativistic limit, for relativistic velocities the friction force fj^2l) is of the form 
Ffj-fA ~ —Vut'jw, which does not saturate in the ultra-relativistic limit. As a consequence, 
this model gives always stationary solutions. Even for i^nr — )■ we have detonation 
solutions and never runaway solutions. 



4.1.2 A friction which saturates 

Recently [13], a simple modification of the usual damping term u^d^j^cj), was considered in 
order to take into account the friction saturation. In our notation, the damping term is 
of the form 

/(0) u^d,cl> 

The constant vector is given by = (0, 0, 0, 1) in the wall frame, and the function 
/(0) is given by Eq. (HTj) . /(0) = (i7nr/T)0^, so that the friction parameter can be 
determined by comparison with microphysics calculations in the non-relativistic regime. 
With the modification the factor '-yv in Eq. (I38p gets replaced by v. This gives a 
friction per unit area of the form 

Fir /A = vMv), (44) 

which has a correct small- velocity behavior as well as Eq. (1421) and, besides, saturates for 
large 7f . More recently, in Ref. |i6j the old model was considered, again with a function 
/(0) of the form /((/>) = (r^/T)^^, but with a velocity-dependent coefficient r] = ?7o/7) so 
that the 7 factor in Eq. fl55]) cancels out. The coefficient r/o is determined by calibrating 
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with non-relativistic results. Therefore, this approach gives again a friction of the form 

dm). 

Although the friction f l44p saturates for f — )■ 1, we note that this approximation is 
too simplistic, as the friction coefficient is determined in the NR limit. Numerically, this 
model cannot give a correct value for the friction in the UR limit, unless t^ur = '7nr- 
Conversely, if the coefficient in Eq. dH]) were determined by ?7ur instead of ?7nr, then the 
model would not give the correct value of the friction in the non-relativistic limit. Notice 
that r/NR and ?7ur will have in general different parametric dependence, given, e.g., by Eqs. 
fll8til9p or by Eq. (l33i) . respectively. It is clear that a model with a single parameter falls 
short of describing the friction in the two opposite regimes. A more realistic interpolating 
model should include two free parameters. 



4.2 Adding a free parameter 

The vector in the model fH3l) may arise as an effective value of a covariant vector. It 
is associated to the presence of the wall and, therefore, should depend on gradients of the 
field. Therefore, we propose the following phenomenological equation for (j), 

,^^,,?r n£iM^ ,45) 

We have chosen a form of the damping term similar to that of Eq. fH5]) . which will 
allow us to compare our results to those of Refs. [13 US]. Notice that other forms [e.g., 
u^d^cp/il + A^M^)] reproduce as well the two behaviors Ffr ~ v and Ffr ~ constant in 
the NR and UR limits, respectively. Thus, in our model the phenomenological vector A^ 
arises from derivatives of the field. In practice, the important difference with the model 
(143 p is that the z component A^ is not set to 1, but gives a new free parameter which can 
be set so as to obtain the correct UR limit. 

Let us consider a reference frame which instantaneously moves with the wall (in which 
most time-derivatives vanish). The phenomenological damping term is given by 



^Jl+[g{ct))(t)'lvX 



The field derivative appearing in the denominator is responsible for a correct correspon- 
dence with the deviations from local equilibrium, as discussed in section [31 In the non- 
relativistic limit, the damping term f HBj) is proportional to cf)' and, thus, is localized inside 
the wall. This corresponds to assuming that the departure from equilibrium occurs inside 
the wall. This is reasonable for a slow wall since, after the wall has passed through a given 
point in space, the system quickly reaches the thermodynamical equilibrium. On the other 
hand, in the ultra-relativistic limit 0' cancels out (together with 7w). For /(0) oc cjP' as in 
Eq. fl4ip . the damping term vanishes in front of the wall but not behind it, since varies 
from 0+ = to 0_ 7^ 0. This is also reasonable since, as we have seen, in the runaway 
regime the deviations from equilibrium occur at the wall and behind. 
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Proceeding as before, in the reference frame which moves with the wall we obtain 



J d(t)dz J ^^TJHWgWWf 

The friction is given by the second integral in Eq. ( H7|) . For f ^ 1 we recover Eqs. 
(138113911. i.e., we obtain a linear friction force. 



A 



= I dzf{(f)){(j)') = -r]t^RV^„ (48) 

NR 



whereas for 7^ ^ 1 we obtain a constant friction force. 



A 



UR 



Md</) = -T^uR. (49) 

9W) 



The functions /((/>) and g{(j)) can be chosen so that Eqs. (HHI) and (H9l) give the desired 
values of the friction coefficients r^NR and t^ur in different models. For instance, a function 
/(0) = r)NR0^/T gives again a friction coefficient of the form (ITH]) in the non-relativistic 
limit. With this choice for f{4>), a 0- independent function 

9W = (50) 

gives a friction coefficient of the form in the ultra-relativistic limit. For intermediate 
velocities the friction depends on the velocity profile as well as on the field profile. 

Solving the set of differential equations for v, T and is out of the scope of the 
present paper. Instead of that, we shall use the thin-wall approximation as before in 
order to obtain analytical results. Invoking again the weighted mean-value theorem for 
integrals, we notice that the friction integral in Eq. f H7|) can be evaluated by replacing the 
whole coefficient of (0')^ with its value at a certain point z = z inside the walo It should 
be a good approximation to assume that (p'i^z) is symmetric and, thus, picks the values 
(p = (f)o/2 and \z = g{(f)o/2)(l)'{z). In the thin- wall approximation one can estimate 0' as, 
e.g., 0' ~ 4>o/lw In fact, whatever approximation we use for 0', we still have the freedom 
to choose the function 5f(0) so as to obtain the desired value of the friction parameter 
r^uR- With these approximations the model essentially reduces to Eq. only that in 

our case Xz 7^ 1. It remains to find an approximation for v{z). The simplest and more 
reasonable one would be v{z) ^ (t>) = (t>+ + t>_)/2. This gives a friction force 

i^fr _ mR{v) ^^^^ 



2 



■^Although the damping term P5|) does not vanish behind the wah in the UR hmit, in Eq. (|T7)) there 
is an extra factor of 0' in the integral for the force acting on the walL 
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Figure 1: A friction force of the form fui/y 1 — (1 — A^)f^. From top to bottom the 
curves correspond to A = 0,0.1,0.2,0.5, 1, 1.5 and 3. Red hues indicate the cases A = 
and A = 1. 

Notice that in Eq. the approximation (^7) was used instead of {v)'y{{v)). In order 
to compare with previous results, in the present paper we shall use, accordingly, 

T-''»(7T=T^) (^2) 

instead of Eq. fl5ip . so that for A^ = we recover Eq. (1421) . For A^ = 1 we recover Eq. 
(144 p . There should not be a significant difference in using either of Eqs. (151 p . (I52p . 

The behavior of this friction force for different values of A^ is shown in Fig. [1] (where 
the effects of hydrodynamics have been neglected, so that v = —Vw). The old model fl37|l 
corresponds to the case A2 = and gives a friction which never saturates. We see that, 
the larger the value of A^, the sooner the friction saturates. The models considered in 
Refs. [TB] correspond to the case A^ = 1. This gives a friction which, in the ultra- 
relativistic limit, is a constant proportional to ?7nr. In the limit v^] — )■ 1, our model gives 
Ff^/A = — r/NR/Az. Comparing with Eq. (l29l) . we see that the parameter A^ is given by 

A. = ^, (53) 

'7UR 

which shows explicitly that the particular case A2 = 1 corresponds to r/uR = V^Ry whereas 
the case A^ = corresponds to the limit t^ur -^00. In terms of the two independent 
friction parameters, the equation for the wall reads 

az^ = p.{T-) - p+(T+) + (.) (T+ - T_) + ( ^====^ ) , (54) 

where we have used the approximation f[T^ for the driving force. 



17 



5 The wall velocity 



In this section we shall study the wall velocity, considering independent variations of the 
friction coefficients, the nucleation temperature T„, and other thermodynamical quan- 
tities. To proceed further, we need to consider an equation of state. It is convenient 
to use the well-known bag EOS, which is simple enough to obtain analytical and model 
independent results. 

5.1 The bag EOS 

The bag EOS can be derived from a free energy density of the form 

jr^ (T) = -a+TV3 + e, J^^ (T) = -a_TV3. (55) 

Thus, the metastable phase is characterized by radiation and false vacuum, whereas in 
the stable phase we only have radiation. In this model the speed of sound is a constant, 
c± = 1/ y/S. More generally, we could consider a nonvanishing vacuum energy density e_ 
in the stable phase. For simplicity we just set e_ = and e+ = e. The only difference is 
that all our results below would otherwise depencQ on Ae = e+ — e_ instead of e. Some 
words are worth, though, on using these results. In a real model, e+ and e_ will be given 
in general by the scale v of the model. Thus, we will naturally have e± ~ f^. In contrast, 
Ae may be much smaller than that, since only a part of the false vacuum energy is in 
general released at T = T^. For instance, in a second-order phase transition we will have 
Ae = 0. Hence, for application of our results, it must be taken into account that Ae = e 
corresponds to the false vacuum energy density that is released at T = Tc, and not the 
total false vacuum energy density. For that reason, it is convenient to use, instead of e, 
the latent heat L, which in a particular model can be easily calculated as explained in 
section |2l Notice that L is larger than e, since thermal energy is released in addition to 
false vacuum energy. For the bag EOS the critical temperature is given by the equation 

(a+ — a-)T^ = 3e, (56) 

and we have the simple relation 

L = 4e, (57) 

since a thermal energy density (a+ — aJ)T^ is released in addition to the vacuum energy 
density e. 

We shall consider bubbles nucleated at a temperature T„ < T^. Hydrodynamics will 
be dominated by the pressure difference between phases, Ap = p^(Tn) — ]9+(T„), and by 
the energy density which is released at the phase transition fronts, Ap = p+(T„) — p_(T„). 
For Tn close to Tc we have Ap ~ and Ap ~ L. In contrast, for T„ — )■ (i.e., for strong 
supercooling), the pressure difference, as well as the released energy, are just given by 

^The false vacuum energy density affects the development of a phase transition through the dependence 
of the expansion rate of the universe on the total energy density. However, our results for the wall velocity 
as a function of thermodynamic and friction parameters depend only on Ae. 
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Ap = Ap = e. As we have seen, in the general case the driving force is not just given by 
Ap, since it is affected by temperature gradients. 

Let us first consider stationary solutions. The fluid discontinuity equations (jUlSD can 
be used to eliminate the pressure difference p_ — p+ from Eq. fl5ll) . For the bag EOS, 
this gives a relatively simple equation |20l 126] . 



2v+V- 



:i + s)(l-f) , r/NR 



1 — 3v+V- 



3a^ 



L 



+ 



where 



a+Tf 4a I Tf ' 



0, (58) 



(59) 



-T, 



and 



a+ / 1 + v+v- ^ 

— 1 - 3aH 

a_ V 1 — 6v+v^ 



1/4 



(60) 



(61) 



In the last expression, we used again the fluid equations (jUlS]) to write the enthalpy ratio 
W-/w+ in terms of v+ and V-. Thus, Eq. fl58|l relates the fluid variables f +, f_, a+. Besides 
these variables, Eq. fl58|) only depends on the parameter ratios //nr/-^, '7ur/''?nr5 
The later can be written as a_/a+ = 1 — Sa^, with Uc = e/{a^T^). Additionally, Eqs. 
(jUlS]) give a relation between and u_ 



1 



1 + an 



6f_ 



± 



1 



6f_ 



Y 



(62) 



The plus and minus signs in Eq. (162 p indicate that we have two hydrodynamical solutions, 
namely, detonations and deflagrations. Using this relation, one can readily solve Eq. fl58|l 
to obtain the velocities v± as functions of a+ (hence, as functions of T+). To obtain the wall 
velocity as a function of the nucleation temperature T„, appropriate boundary conditions 
must be used. The relation between the fluid velocities v± and the wall velocity f^, as 
well as the relation between a;+ and a;„ = e/ (a+T^), depend on the type of hydrodynamic 
solution. For detonations we have simply a+ = and f+ = —Vw For deflagrations, the 
matching conditions at the shock discontinuity must be used. These give 



Vw - \V+\ 

1 - \v+\v,, 



a/3 (a„ - a+) 



a/ (3a„ + a+) (3a+ + 



(63) 



whereas for Jouguet deflagrations we 



For traditional deflagrations we have f_ = — 
have f _ = Co = — l/VS (for details see Refs. [2D1ES]). 

Using the bag EOS, we can also express the runaway conditions in terms of the ther- 
modynamic parameters and the ultra-relativistic friction coefficient. From Eqs. f l55ll5T 
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we have p^(Tn) — p+(Tn) = (L/4)(l — T^/T^), and the necessary condition (l35i) gives 
simply 

^ + — <1- (64) 
We may express this condition as t^ur < Vnec, with 

^nec/e = 1 - ac/a„. (65) 
In terms of the a variables, we have 

an > anec = "j — — — T" (66) 

We remark that Eqs. ( I64ll66l) only give a necessary condition for the wall to run away, as 
the runaway solution may coexist with a detonation solution, which is presumably stable. 
Nevertheless, detonations are not possible if the sufficient condition ( !36|) is fulfilled. For 
the bag EOS and for a detonation with ~ 1, we have p_(T_) — p_|_(T+) = L/2, T_/T^ = 
(a_/a+)~^/^(l + 3a„y/^, and s_/s+ = {a_/a+Y^'^{l + 3anY^^ (for analytic approximations 
for ultra- relativistic detonations see Ref. [26]). We thus obtain the sufficient runaway 
condition r/uR < Vsui, with 



Vsuf 



3ar. 



;i - 3a,)3(l + 3ajt - (1 - 3a,)-3(l + 3a„)i 



(67) 



Equivalently, for given values of rjuji/e and etc we may express this condition as > Osuf, 
i.e., in terms of the amount of supercooling which is sufficient for the wall to runaway. 
The value of Oguf is obtained by inverting Eq. (jHT]). This amounts to solving a quartic 
equation for the variable x = (1 + 3a„)^/^. The expression for Oguf is cumbersome and we 
shall not write it down. 

Before going on to the interpretation of these results, it is worth comparing them 
with the previous works [TB] . The runaway necessary condition can be written in 
terms of fundamental parameters rather than phenomenological ones. Going back to the 
original expression p_(T+) > p_|_(T+), and taking into account the fact that p_ is given 
by the 0((/>^) expansion of the thermal part of the effective potential, one may use the 
approximation fl3T|) to write the condition as e > T^0?^ ^Cj(7j/i^/24 or, dividing by a+T^, 
as 

Of course, this condition is equivalent to Eq. fIBB]) . i.e., a^ec = ttc/(l — Vur/^)- Although 
Eq. (168|) is discussed in Ref. [15], the simple phenomenological model used for explicit 
calculations corresponds to setting ?7ur = tj^r, which gives a^ec = ac/(l — V^r/^)- Sim- 
ilarly, in Ref. [I6j the dependence of the microphysics on the wall velocity is discussed 
in some detail. However, when hydrodynamics is included in the calculation, a friction of 
the form rjo/'y is assumed. As we have already mentioned, this phenomenological model 
corresponds again to r^uR = '7nr- 
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Figure 2: Regions in the (T„, ?7uR)-plane where runaway and detonation solutions can 
exist. The runaway necessary condition f lM|l is fulfilled below the blue lines. The sufficient 
condition (|67|) is fulfilled below the red lines. 



5.2 Stationary and runaway regimes 

From Eq. f l64|) we see that the existence of runaway solutions requires that both the UR 
friction parameter and the nucleation temperature be small enough. In particular, the 
necessary condition is never fulfilled for r/uR > ^i no matter how small Tn/T^.. This is 
because the pressure difference Ap(T„) is bounded by e (Ap — )■ e for T„ — )■ 0) as much as 
the friction is bounded by t^ur- For ^ur smaller than, but close to e, a strong supercooling 
is required [i.e., (T„/Tc)^ <C 1]. Conversely, if the amount of supercooling is small (i.e., 
Tn ~ Tc), then ?7ur <^ e is required. The quantity e = L/4 is related to the strength 
of the phase transition. Notice that the parameter = L/(4a+T^^) is limited to the 
range < < 1/3. The lower limit corresponds to a second-order phase transition with 
a_ = a+ and L = 0. The higher limit corresponds to the case a_ = and L = TcS+(Tc), 
i.e., to maximum entropy discontinuity at T = Tc. Figure |2] illustrates the necessary 
condition fl64|) and the sufficient condition f l67|) in parameter space. In terms of rjim/e, 
the necessary condition does not depend on the parameter ac (left panel). In terms of 
rju^/{a+T^), it does (right panel). 

In the region above the blue lines the runaway solution does not exist. Below these 
curves the runaway solution exists, but it may coexist with stationary solutions. The red 
lines indicate the limits of existence of detonation solutions. Thus, for a given value of ac, 
the region above the corresponding red curve belongs to detonations or other stationary 
solutions, whereas the runaway solution is possible, in principle, only below the curve. 
Qualitatively, the two conditions have a similar behavior in the (T„, 77uR)-plane. Both 
require small enough values of t^ur and T„. Quantitatively, the sufficient condition is 
more restrictive and has a stronger dependence on hydrodynamics. This is because, for 
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the stationary solution, the released energy causes reheating and bulk motions of the fluid. 
The runaway wall causes less perturbations, since the released energy goes mainly into 
accelerating the wall. 

Runaway solutions are driven by the pressure difference Ap(T„). For close to Tc, 
we have Ap(T„) ^ and there won't be runaway solutions unless the friction parameter 
is very small, as can be seen in Fig. [2l As we increase the amount of supercooling, the 
pressure difference increases. For a given value of rju^i/e, one may in principle expect that, 
as we keep increasing the amount of supercooling, runaway solutions will eventually be 
possible. As we have seen, though, the pressure difference saturates for strong supercooling 
(as does the friction force for high wall velocity). Thus, the wall will not run away for 
^UR > e (see the left panel). Below the blue curve, runaway solutions become possible. 
However, whether the wall will run away or not, depends on the values of the parameters. 
Consider for instance the case etc = 0.1 (dotted hne) and a fixed value of t^ur- If the 
value of rju^/e is, say, 0.8, then the wall will never run away, no matter how strong the 
supercooling. For a smaller friction, say, //ur/^ = 0.6, then the wall will run away for a 
strong enough supercooling (T„/Tc < 0.66). Eventually, for a very strong supercooling 
(Tn/Tc < 0.29) we recover again the stationary solution. 

This behavior is due to the added effects of friction and hydrodynamics. Consider a 
weakly first-order phase transition, which is characterized by a small amount of released 
energy (i.e., L <^ a+T^, which implies ac <^ 1) and a small amount of supercooling (i.e., 
Tn ~ Tc). In such a case, the effects of hydrodynamics will disappear and the sufficient 
condition will approach the necessary condition, as can be clearly seen in Fig. [2] (see, e.g., 
the cases Oc = 10~^ and ac = 10~^). A strong phase transition is characterized by a large 
latent heat as well as a significant amount of supercooling. For T„ ^ Tc (and t^ur < 
one expects that the wall will run away. However, the release of latent heat slows down the 
detonation, since the driving force is affected by temperature gradients. For a fixed value 
of the ratio L/ (a+T^) this effect becomes more and more important as the released energy 
density Ap ~ T^^ becomes larg^ in comparison with the plasma energy density p ~ T^. 
This is why the sufficient condition departs further from the necessary condition for small 
Tn. The left panel of figure [2] shows that, for extremely supercooled phase transitions, the 
friction must be very small for the wall to run away. 

In fact, microphysics gives temperature-dependent friction parameters. In a realistic 
model, we expect that ?7ur will decrease at small temperatures, since the density of par- 
ticles in front of the wall vanishes for T — )■ 0. Hence, the right panel in Fig. [2] may be 
more illustrative for small temperatures. Although the curves are more involved (since the 
necessary condition now depends on ac) we see that, if we fix the value of 77ur/('3-+T^), 
the system will always enter the runaway region for strong enough supercooling. We 
shall discuss the strong supercooling case (as well as the actual dependence of friction on 
temperature) in Sec. O 

The behavior of the red curves can also be seen analytically. For a„ ~ etc <^ 1, 
Eq. ([HI]) gives ?7suf/e = r/ncc/e - (3/8)(ac/an)(an + "c) + 0{al). The negative sign of 

^For the bag model, the released energy density is bounded below by e = L/4. 
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the 0{a) terms implies that the value rjsui is smaller than rj^ec- As we increase etc, the 
first term in Eq. (I^7|) decreases whereas the second term increases. Hence, the value of 
r^suf becomes smaller, departing from T^nec- As a consequence, the region where the wall 
can runaway gets reduced. In the strong supercooling limit, we have a linear behavior, 
r/suf/e = (2/3i/^)(lK - Sf/^TJ^) + 0(T„/T,)3. 

5.3 Stationary wall velocity 

We shall now study the stationary wall velocity as a function of the parameters. We shall 
start the parameter variation from previously considered values [121 [131 [HI [20], which 
are convenient to exhibit all the kinds of hydrodynamical solutions and to compare with 
previous results for the old phenomenological friction model. In Fig. [3] we considered 
the stationary wall velocity as a function of the non-relativistic friction parameter t^nr- 
Although we can vary independently ?7nr and r/uR? in a real model there will be some 
correlation between these parameters (e.g., both increase with the couplings of particles 
to the Higgs). Therefore, we considered two opposite relations for them. In the upper 
panels of Fig. [3] we fixed the ratio = rj^R/rjuR, whereas in the lower panels we fixed 
?7uR- The left panels show all the stationary solutions for a given set of parameters. The 
right panel shows the solutions which are actually realized in the phase transition. 

Consider first the case of fixed (upper panels). In the left panel we see that, in some 
cases, there are several possible solutions. In the first place, detonations and "traditional" 
deflagrations are bivalued in a certain range. This is not a problem, since the lower 
branch of detonations and the upper branch of deflagrations are unphysical solutions 
which, according to numerical calculations, are unstable [121 [131 [Il| and, thus, must 
bTdiscardecfi. On tbe other hand, we see that also different types of hydrodynamical 
solutions may coexist. Numerical calculations seem to indicate that supersonic traditional 
deflagrations, which are strong deflagrations, are unstable. Thus, the only supersonic 
deflagration which can be realized is the Jouguet deflagration. In case of coexistence of a 
deflagration with a detonation, the detonation seems to be the stable solution [121 [131 [SI • 
We have plotted the curves in the right panel assuming this hierarchy. 

The case A^ = (solid line) corresponds to the widely used model (137|) . which does 
not exhibit runaway behavior. For Xz ^ the friction (l52i) saturates at high velocities 
and the wall may run away. The runaway behavior depends only on ijur. However, 
for fixed = VNR/vuRy the wall always runs away for small enough ?7nr, i.e., the wall 
velocity becomes f ^ = 1 at ?7nr = ?7suf / Xz (indicating that the stationary wall assumption 
breaks down). Notice that the behavior is qualitatively similar for the different values 
of A^. Quantitatively, the differences become significant only for small ?7nr or large Xz- 
The particular case A^ = 1 corresponds to the models of Refs. [151 [IB], for which the 
friction coefficient has the same value in the two opposite regimes. Lower values of A2, 
corresponding to higher values of the UR friction, yield a wider range of parameters with 
stationary solutions. 

^Notice that this ehminates the Jouguet detonation (corresponding to the lower endpoint of the 
detonation curves), which has often been considered for gravitational wave generation. 
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Figure 3: The wall velocity for the bag model with ac = 4.45 x 10"^ and T„ = 0.89Tc. 
Detonations are in blue, "traditional" deflagrations are in red, and Jouguet deflagrations 
are in black. The right panels show only the stable solutions. 
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Figure 4: The wall velocity as a function of ?7ur, for the same bag parameters of Fig. |3] 
and several values of t^nr- The vertical lines indicate the values of r/guf and rjnec- 



In the lower panels of Fig. |3] we considered fixed values of t^ur- We have chosen 
most values around r^uR = '7suf — 0.0904 (dashed-dotted line). For higher values of r^uR 
the curves accumulate near the limiting line of r/uR ~^ oo (which is the same as the 
curve of A2 = in the upper figure). For t^ur < Vsuf the detonation solutions disappear 
(the curves in the left panel correspond to the unphysical branches of the solutions, for 
which the velocity increases with the friction), and the deflagrations approach the speed 
of sound. The right panel shows only the physical stationary solutions. It is interesting 
that for r^uR < ''7suf there are still deflagration solutions. We show these curves in grey. 
However, these fast deflagrations coexist with the runaway solution, and it is possible 
that they are unstabl^. Consequently, we shall assume that for r^uR < Vsui the wall runs 
away. For r/uR > Vsuf, we notice that there is always a physical stationary solution, no 
matter how small the value of the NR friction coefficient. This behavior contrasts with 
the upper panel, where the wall always runs away for small enough t^nr. As we have 
already mentioned, in a physical model there will always be a correlation between i^nr 
and ?7uR. We do not expect, in general, a very small t^nr together with a considerably 
large ?7ur or vice versa. As the fundamental parameters of a model are varied, we expect 
an intermediate behavior between the curves of the upper and lower panels. 

In Fig. IHwe flxed ?7nr and varied t^ur- Notice that the wall velocity changes quickly 
near the critical value ?7ur = ''^suf- For larger values the velocity is essentially a constant 
which depends on the value of ?7nr. As t^ur approaches ?7suf from the right, the velocity 

''in general, the traditional deflagration seems to become unstable when it coexists with other station- 
ary solutions (namely, Jouguet deflagrations or detonations) . It is probable that the same happens when 
it coexists with the runaway solution. A stability analysis is out of the scope of this paper and shall be 
addressed elsewhere. 
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Figure 5: The wall velocity as a function of an/o-c = {Tc/TnY, for etc = 4.45 x 10 ^, 
titssr/ L = 0.1, and several values of A^. 

grows and the detonation solutions disappear (their velocity becomes v^, = 1). Again, 
we see that for r/uR < Vsuf we may still have physical deflagration solutions (in gray in 
the right panel). As explained above, we shall assume that these are unstable and we 
shall chose the runaway solution. In the right panel we have zoomed the friction near 
^UR = Vsui and we have marked the values of ?7suf and rj^ec (which are very close to each 
other for the present values of ac and 

Finally, in Fig. Elwe fixed the friction parameters and varied the amount of supercool- 
ing. In the left panel we also show the speed of sound and the Jouguet detonation velocity 
f j^*(a„) (dotted grey lines). We see that, for T„ close to Tc, the solution is always a weak 
deflagration and does not depend on the UR friction parameter. The detonation solu- 
tions reach the speed of light at a„ = Osuf ? indicating that beyond that value the wall runs 
away. The runaway solutions exist already from a„ = anec, indicated in the right panel 
by black round dots on the curves. For the smallest values of (highest values of //ur) 
the velocity grows very slowly with the amount of supercooling. As we have discussed 
earlier, the wall may eventually decrease for strong supercooling, due to hydrodynamics 
effects. We shall analyze this possibility in the next section. For high values of A^, the 
ranges of stationary solutions get reduced (see the right panel). This is also observed in 
the previous figures. For the case A^ = 2, the Jouguet deflagration has disappeared (the 
weak deflagration turns directly into a detonation) and the detonation has a very short 
parameter range of existence. For Az = 5 the detonation has also disappeared and the 
weak deflagration turns directly into a runaway solution at a„ = Osui (we have marked 
this value with a squared blue dot, and we have plotted the deflagration solution in grey 
beyond it). As can be seen in the left panel, in this case it is the unphysical branch of 
detonations that reaches the speed of light. 
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In all these figures, we see that the detonation velocity becomes w^^, = 1 at the values 
Osuf or ?7suf obtained in sec. 15.11 As we have just seen, in some cases this actually 
corresponds to an unphysical detonation. In general, cases like that of = 5 in Fig. 
[5l where the physical detonation does not exist at all, deflagrations do exist. Assuming, 
as we did, that in this case the deflagration becomes unstable and jumps to a runaway 
solution at a„ = Osuf (or t^ur = Vsuf in the right panel of Fig. H]), the values of ctguf 
and risui still make sense as limits between a stationary solution (either a detonation or a 
deflagration) to a runaway solution. This is the actual meaning of the red curves in Fig. 

El 

6 Strong supercooling 

The fast propagating modes of a phase transition front, namely, detonations and runaway 
walls, are interesting due to the possibility of generating sizeable gravitational waves. The 
high velocity behavior is governed by the UR friction parameter. As we have seen, though, 
whether the wall will be fast or not, depends also on the amount of supercooling and on 
the latent heat. A considerable amount of supercooling favors large pressure differences 
and high velocities, whereas a large release of latent heat slows down the wall. Strong 
supercooling is typical of strongly first-order phase transitions. However, the latter are 
also characterized by a large latent heat. On the other hand, at low temperatures one 
expects a small friction. 

The friction depends strongly on the values of the particle masses inside the bubble 
(assuming for simplicity that rrii = outside); more precisely, on the ratio rrii/T. Thus, 
a relevant parameter is the value of 0/T in the broken-symmetry phase. Strongly first- 
order phase transitions are characterized by a relatively high value of the ratio (j)o{Tc)/Tc. 
This implies also a relatively high latent heat. Nevertheless, 0o and Tc are both given 
by the characteristic scale of the theory, and its ratio will be 0{1) for a natural phase 
transition. As we have seen, for the bag EOS, L is bounded by ^a+T^. On the other hand, 
in some models the barrier between minima may persist at low temperatures, causing a 
large amount of supercooling. In such a case we may have <^ Tc, which implies 
MTn)/Tn > 1 and L > a+T^. 

On the whole, for very strong supercooling we will have in principle a high pressure 
difference between phases and a small friction force, but the energy released will be large 
in comparison to the energy density of the plasma, and hydrodynamics effects will be 
important. It is thus worth investigating the behavior of phase transition fronts in such 
a case. 

6.1 Detonations and runaway solutions 

As we have seen, strong hydrodynamics effects may prevent the wall from reaching ultra- 
relativistic velocities, even in the case of a strong supercooling. To illustrate this effect, 
let us plot again vs. a„ like in Fig. |5l this time for a wider range of the supercoohng 
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Figure 6: The wall velocity as a function of O-njoLc for rijL fixed (solid) and ri/{aj^T^) 
fixed (dashed). In the left panel we have t^nr = O.IL for T„ close to T^. In the right panel 
we have t^nr = 0.15L and r/uR = 0.175L for close to T^. 

parameter an/oLc- It is interesting to consider first the same thermodynamic parameters 
of Fig. [5l which correspond to a relatively low value of the released energy (since L ~ 
lO'^a+T^"^). In the left panel of Fig. [6] we consider a couple of representative values 
of \z- As an example of the case ?7ur > e, for which the wall will not run away, we 
choose \z = 0, i.e., the limit r^uR oo (black lines). We see that, for fixed values 
of the friction parameters (solid black line), the velocity grows slowly with the amount 
of supercooling, and eventually decreases for very strong supercooling. Moreover, for 
anjoLc > 300 the stationary solution cannot exist due to strong hydrodynamics effects. 
Before disappearing, the detonation approaches the Jouguet point (dotted grey line) and 
a non-physical branch of solutions appears, just like at the other end of the curve (where 
the supercooling becomes insufficient for a detonation solution). As we have mentioned, 
for strong supercooling it is more realisfJl to consider a friction which decreases as T„ 
decreases. Therefore, we have also considered fixed values of ?7/(a+T^) (dashed black 
line). In this case, the wall velocity is only bounded by the relativistic limit = 1- The 
wall does not runaway because we are considering the case ?7ur oo. We also consider 
the case = 0.7 (blue lines). The solid line {rj/L fixed) corresponds to the dotted line 
in Fig. [5l which runs away for anjoic — 2.4. For ri/{a+T^) fixed (dashed blue line) the 
wall runs away at a smaller a„. 

Hydrodynamics effects will be stronger for higher values of the latent heat. In Sec. 15.21 
we have discussed the case ac = 0.1, corresponding to the dotted line in Fig. [21 According 
to that figure (see the left panel), for rju^/L ^ 0.15 (i.e., //ur/^ ~ 0.6) the detonation 
solution gives way to the runaway solution at a given amount of supercooling. But then, 

^We discuss this issue below. 
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for a stronger amount of supercooling, the detonation solution is possible again. The right 
panel of Fig. shows this effect (we have only kept the physical solutions and we have 
chosen the values rj^ji/L = 0.15 and rjuji/ L = 0.175, which better illustrate the effect). 
It is important to notice, however, that shortly after reappearing, the detonation ceases 
to exist due to strong hydrodynamics. In many cases, the detonation will not reappear 
at all. In any case, if we consider a friction which vanishes at zero temperature (dashed 
line), then a smaller supercooling will suffice for the wall to run away, and the stationary 
solution will not appear again at stronger supercooling. This can be seen already in the 
right panel of Fig. |2j 

So far we have considered two simple cases for the variation of the friction coefficients 
with temperature, which consist of fixing the dimensionless ratios t]/L and ri/{a+T^). 
In a realistic model, we expect temperature-dependent friction parameters, although not 
necessarily decreasing as fast as T^. Roughly, the number density of massless particles in 
the symmetric phase goes as T^, and the force each particle feels at the wall is proportional 
to Ap ^ (m/p)Am (in the WKB approximation). For Am = m(0o), we have Ap ~ vn? jT . 
Thus we obtain a force of the form F ~ rp?T'^ . As we shall see, this is parametrically 
correct for r/uR- The case of //nr is more complex because one must take into account 
particle interactions. For interaction rates F ~ T, one obtains again Ffj. oc (see below). 

To investigate the behavior of the wall velocity with a friction which depends quadrat- 
ically on temperature, we shall consider friction parameters of the form ri = rjoT^ = 
vIq^/ ad OLn- We choose the coefficient r^o so that the different forms of the friction match 
for T„ ~ Tf.. In order to evaluate to what extent fast solutions depend on the NR fric- 
tion parameter, we shall also consider the case of a constant r^NR with a quadratic t^ur- 
We show the result in Fig. [TJ together with the curves of Fig. E] for comparison. The 
black lines correspond to both friction coefficients of the form ri = rjoT'^. The red lines 
correspond to a constant NR friction coefficient. We see that the black and red curves 
are very close to each other, and meet in the limit = 1. This is because the value 
of the NR friction coefficient does not play a relevant role for detonations, and becomes 
irrelevant in the ultra- relativistic limit. In the right panel the effects of hydrodynamics 
are stronger since the released energy is larger (notice the larger amounts of supercooling 
needed to obtain detonations and runaway walls). We see that the solid line deviates 
significantly from the others. This corresponds to the unrealistic case in which the UR 
friction is temperature-independent. 

6.2 Strong supercooling and friction coefficients 

In this last subsection we wish to discuss briefly the behavior of the friction coefficients at 
low nucleation temperatures (see |28j for other limiting cases). A detailed analysis is out of 
the scope of the present paper. Therefore, we shall assume that the WKB approximation 
is still valid. For very low temperatures, such that T ^ l^^ this approximation will 
certainly break down. 

Let us first consider the non-relativist ic friction coefficient r^NR- Since the interaction 
rates are proportional to temperature, at low temperatures the approximation F ^ 
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which leads to Eq. (|T6|) is no longer valid. Nevertheless, we shall consider the approxi- 
mation (|T6|) to catch a glimpse of the behavior of the friction. In any case, we have just 
seen that the wall velocity depends very weakly on t^nr for strong supercooling. In the 
limit m/T ^ 1 the coefficient Ci is suppressed by a Boltzmann factor. 



ci ^ (m/T)^/^exp (-m/T) / (27r) 



3/2 



(69) 



However, in the case of interest, very light particles in the symmetric phase become 
heavy in the broken-symmetry phase. For simplicity, let us consider masses of the form 
mi{(j)) = hi(j). To calculate the integral in Eq. f lT6|) . we notice that for small the 
integrand is suppressed by powers of cf). Therefore we can use the approximation fl69l) in 
all the range of integration. We have 



??NR 



r(27r)= 



dx 



dz 



dx, 



(70) 



with X = hi(j)/T. To perform the integral, an approximation for (f){z) is needed. Instead of 
the usual (0o/2)([l + tanh(z//^„)], we shall use the much simpler approximation of a linear 
function (j) = (p^/l^ z inside the wall and = constant outside the wall. The temperature 
dependence, as well as the parametric behavior, should not depend on the approximation 
for the wall profile. Thus, in the integral ( I70|) dx/dz is a constant and we obtain, for 
Ho/T > 1, 

E^gihiCpoT"^ 
^ 64vr3(r/T)L- 

We see that the friction coefficient vanishes for small T, though not as fast as T^. 
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The ultra-relativistic friction coefficient is very easy to estimate in the strong super- 
coohng hmit, either from its definition in Eqs f l28tl2^ or directly from the total force f l20l) . 
It is interesting to consider the latter. For 0+ = and rrii = hi(j), the total force is given, 
at any temperature, by 



At high temperatures, the last term in Eq. fl72]) is part of the equilibrium pressure 
difference [as can be seen in the high-temperature expansion (l3Ti) ]. Hence, it is part of the 
driving force. As a consequence, the friction is given by the rest of the high-T expansion, 
VvR ~ ^0o + ^(0o) ^'i- dS])]- At very low temperatures, on the contrary, the driving 
force is just given by the vacuum potential, i.e., the ffist two terms in Eq. ( 1721) . Hence, 
the last term now gives the friction force, and we have 



We see that r/uR goes as at low temperatures. 

7 Conclusions 

In this paper we have studied the friction force acting on phase transition fronts. In 
particular, we have discussed the ultra-relativistic behavior of the friction, derived in Ref. 
[To] . We have considered runaway walls as well as stationary walls (detonations), which 
have different hydrodynamics. We have shown that the two solutions coexist in a range of 
parameters. Thus, we have argued that the runaway condition given in Ref. [10] is only 
a necessary condition and does not guarantee that the wall will run away. We have found 
a sufficient condition for the wall to run away by using the criterion of the detonation 
reaching the speed of light. 

We have also proposed a phenomenological model for the friction, which interpolates 
between the non-relativistic and the ultra-relativistic behaviors. The main improvement 
of this model with respect to previous ones is the incorporation of a new free parameter 
governing the ultra-relativistic limit of the friction force. The value of this parameter 
can be easily calculated for any specific model. Therefore, our model not only gives a 
saturating friction, but also gives the correct quantitative behavior of the friction in the 
UR limit. This phenomenological model is consistent with the necessary and sufficient 
runaway conditions, which do not depend on the non-relativistic friction coefficient. 

Using the bag equation of state and our phenomenological model for the friction, we 
have studied the runaway conditions, as well as the stationary regime, as functions of the 
thermodynamic and friction parameters. Regarding the runaway conditions, the general 
result is that they are quantitatively similar for small values of the latent heat but depart 
otherwise. If we, for instance, increase the amount of supercooling (for a given value of 




(72) 




(73) 
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the friction), then, in the case of small latent heat, the sufficient condition is reached 
shortly after the necessary condition. On the contrary, for larger latent heat, the larger 
release of energy slows down the phase transition fronts. This causes the detonation to 
persist for larger amounts of supercooling, even if the runaway solution already exists. 
Regarding the stationary wall velocity, we have found, as expected, that for parameters 
which favor a small wall velocity (i.e., for low supercooling, large latent heat, or high 
non-relativistic friction) the solution depends very weakly on the UR friction coefficient. 
Conversely, for parameters which favor fast stationary solutions, the velocity does not 
depend significantly on the NR friction coefficient. 

We have studied in particular the case of a phase transition with a large amount of 
supercooling, which is important for the generation of gravitational waves. In this case the 
dependence of the friction on temperature becomes relevant. We have explored the wall 
velocity for different behaviors of the friction coefficients. Thus, we have considered fixed 
values of the dimensionless parameters r]/L and rj/T^, as well as the quadratic dependence 
7] cc T"^. The latter is motivated by low-temperature approximations, which we have also 
discussed. The value of the NR friction parameter does not affect the strong supercooling 
case, since the velocity is in general high. The general result is that, if the UR friction 
parameter does not decrease for strong supercooling, then the wall may reach a saturation 
value, and even decrease due to hydrodynamics effects. In contrast, if the UR friction 
parameter vanishes at zero temperature (which seems to be the actual case), then the 
wall always runs away for strong enough supercooling. 

Acknowledgements 

This work was supported by Universidad Nacional de Mar del Plata, Argentina, grant 
EXA 607/12. 

References 

[1] For reviews, see A. G. Cohen, D. B. Kaplan and A. E. Nelson, Ann. Rev. Nucl. 
Part. Sci. 43, 27 (1993) |arXiv:hep-ph793022iO|; A. Riotto and M. Trodden, 
Ann. Rev. Nucl. Part. Sci. 49, 35 (1999) |arXiv:h ep-ph/9901362] ; T. Konstandin, 
larXiv: 1302.6713 [hep-ph]. 

[2] For a review, see D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001) 
[arXiv:astro-ph /0009061|. 

[3] A. Vilenkin and E.P.S. Shellard, Cosmic Strings and Other Topological Defects (Cam- 
bridge University Press, Cambridge, England, 1994); A. Vilenkin, Phys. Rept. 121, 
263 (1985). 

[4] E. Witten, Phys. Rev. D 30, 272 (1984); G. M. Fuller, G. J. Mathews and C. R. Al- 
cock, Phys. Rev. D 37, 1380 (1988); J. H. Applegate and C. J. Hogan, Phys. 



32 



Rev. D 31, 3037 (1985); H. Kurki-Suonio, Phys. Rev. D 37, 2104 (1988); J. Ig- 
natius, K. Kajantie, H. Kurki-Suonio and M. Laine, Phys. Rev. D 50, 3738 (1994) 
I arXiv: hep-ph/9405336] ; A. Megevand and F. Astorga, Phys. Rev. D 71, 023502 
(2005). 

[5] A. F. Heckler, Phys. Rev. D 51 405, (1995) tarXiv:astro-ph/9407064| ; 

[6] A. Kosowsky and M. S. Turner, Phys. Rev. D 47, 4372 (1993); M. Kamionkowski, 
A. Kosowsky and M. S. Turner, Phys. Rev. D 49, 2837 (1994); A. Kosowsky, A. Mack 
and T. Kahniashvih, Phys. Rev. D 66, 024030 (2002); A. D. Dolgov, D. Grasso and 
A. Nicohs, Phys. Rev. D 66, 103505 (2002); C. Caprini and R. Durrer, Phys. Rev. 
D 74, 063521 (2006); C. Caprini, R. Durrer and G. Servant, Phys. Rev. D 77, 
124015 (2008) j^rXiv:071 1 . 25931 [astro-ph] ] : R. Apreda, M. Maggiore, A. Nicohs and 
A. Riotto, NucL Phys. B 631, 342 (2002); A. Nicohs, Class. Quant. Grav. 21, L27 
(2004); C. Grojean and G. Servant, Phys. Rev. D 75, 043507 (2007); S. J. Huber and 
T. Konstandin, JCAP 0809, 022 (2008) [arXiv:0806.1828 [hep-ph]]; S. J. Huber and 
T. Konstandin, JCAP 0805, 017 (2008) [arXiv:0709.2091 [hep-ph]]. A. Megevand, 
Phys. Rev. D 78 (2008) 084003 [arXiv:0804.0391 [astro-ph]]; J. Kehayias and S. Pro- 
fumo, JCAP 1003, 003 (2010) [arXiv:0911.0687 [hep-ph]]; L. Leitao, A. Megevand 
and A. D. Sanchez, JCAP 1210, 024 (2012) [arXiv:1205.3070i [astro-ph. CO]]. 

[7] B. H. Liu, L. D. McLerran and N. Turok, Phys. Rev. D 46, 2668 (1992); N. Turok, 
Phys. Rev. Lett. 68, 1803 (1992); M. Dine, R. G. Leigh, P. Y. Huet, A. D. Linde 
and D. A. Linde, Phys. Rev. D 46, 550 (1992) fi?Xrv:hep-ph/9203203|; S. Y. Khleb- 
nikov, Phys. R ev. D 46, 3223 (1992); P. Arnold, Phys. Rev. D 48, 1539 (1993) 
arXiv:hep-p h/9302258 1 ; G. D. Moore and T. Prokopec, Phys. Rev. D 52, 7182 (1995) 
arXiv:hep- ph/9506475|; Phys. Rev. Lett. 75, 777 (1995) |arXiv:hep- ph/9503296|; 
G. D. Moore and N. Turok, Phys. Rev. D 55, 6538 (1997) [ar Xiv:hep-p h/9608350|. 

[8] P. John and M. G. Schmidt, Nucl. Phys. B 598, 291 (2001) [Erratum-ibid. B 648, 
449 (2003)]. 

[9] G. D. Moore, JHEP 0003, 006 (2000); 

[10] D. Bodeker and G. D. Moore, JCAP 0905, 009 (2009) [arXiv:0903.409 9' [hep-ph]]. 

[11] M. Gyulassy, K. Kajantie, H. Kurki-Suonio and L. D. McLerran, Nucl. Phys. B 237, 
477 (1984); K. Enqvist, J. Ignatius, K. Kajantie and K. Rummukainen, Phys. Rev. 
D 45, 3415 (1992); H. Kurki-Suonio, Nucl. Phys. B 255, 231 (1985); K. Kajantie 
and H. Kurki-Suonio, Phys. Rev. D 34, 1719 (1986); 

[12] J. Ignatius, K. Kajantie, H. Kurki-Suonio and M. Laine, Phys. Rev. D 49, 3854 
(1994); 

[13] H. Kurki-Suonio and M. Laine, Phys. Rev. D 51, 5431 (1995) 
[arXiv: hep-ph/950 1216] ; 



33 



[14] H. Kurki-Suonio and M. Laine, Phys. Rev. D 54, 7163 (1996) |hep-ph/ 9512"202]. 

[15] J. R. Espinosa, T. Konstandin, J. M. No and G. Servant, JCAP 1006, 028 (2010) 
[arXiv: 1004.41871 [hep-ph]]; 

[16] S. J. Huber and M. Sopena. [a?XwIT302.1044 [hep-ph]. 



[17] M. Quiros, arXiv:hep-ph/9901312 



[18] S. R. Coleman, Phys. Rev. D 15, 2929 (1977) [Erratum-ibid. D 16, 1248 (1977)]; 
C. G. Callan and S. R. Coleman, Phys. Rev. D 16, 1762 (1977); I. Affleck, Phys. 
Rev. Lett. 46, 388 (1981); A. D. Linde, Nucl. Phys. B 216, 421 (1983) [Erratum-ibid. 
B 223, 544 (1983)]; Phys. Lett. B 100, 37 (1981). 

[19] A. H. Guth and E. J. Weinberg, Phys. Rev. D 23, 876 (1981); G. W. Anderson 
and L. J. Hall, Phys. Rev. D 45, 2685 (1992); A. Megevand, Int. J. Mod. Phys. 
D 9, 7 33 (2000j_^[liep-ph/0006177|; A. Megevand, Phys. Rev. D 64, 027303 (2001) 
; A. Megevand and A. D. Sanchez, Phys. Rev. D 77, 063519 (2008) 
[hep-ph]]. 



hep-ph/0011019 



arXiv:0712.1031 



[20] A. Megevand and A. D. Sanchez, Nucl. Phys. B 865, 217 (2012) jarXiv:1206.2339. 
[astro-ph.CO]]. 

[21] L. Leitao and A. Megevand, Nucl. Phys. B 844, 450 (2011) |arXiv:1010.2"T34l [astro- 
ph.CO]]. 

[22] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Pergamon Press, New York, 1989); 
R. Courant and K. O. Friedrichs, Supersonic Flow and Shock Waves (Springer- Verlag, 
Berlin, 1985). 

[23] P. J. Steinhardt, Phys. Rev. D 25, 2074 (1982). 

[24] B. Link, Phys. Rev. Lett. 68, 2425 (1992); K. Kajantie, Phys. Lett. B 285, 
331 (1992); M. Kamionkowski and K. Freese, Phys. Rev. Lett. 69, 2743 (1992) 
[arXiv:hep-ph/9208202|; P. Y. Huet, K. Kajantie, R. G. Leigh, B. H. Liu and 
L. D. McLerran, Phys. Rev. D 48, 2477 (1993) | arXiv:hep-ph/92 12224] . 

[25] A. Megevand, Phys. Rev. D 69, 103521 (2004). 

[26] A. Megevand and A. D. Sanchez, Nucl. Phys. B 820, 47 (2009) iarXiv:0904. 17531 
[hep-ph]]. 

[27] T. Konstandin and J. M. No, JCAP 1102, 008 (2011) [ arXiv: 101 1.37351 [hep-ph]]. 

[28] A. Megevand and A. D. Sanchez, Nucl. Phys. B 825, 151 (2010) |arXiv:0908.3"663l 
[hep-ph]]. 



34 



[29] S. J. Huber and M. Sopena, Phys. Rev. D 85, 103507 (2012) |arXiv:1112.1"888l 
ph]]. 

[30] T. M. Apostol, Calculus, Vol. I, (John Wiley & Sons, New York, 1967). 



35 



